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SUMMARY 


A Damage-threshold/Fai 1-saf ety approach is proposed for ensuring that 
composite structures are both sufficiently durable for economy of operation, as 
well as adequately fail-safe or damage tolerant for flight safety. Matrix cracks 
are assumed to exist throughout the off-axis plies. Delamination onset is 
predicted using a strain energy release rate characterization. Delamination 
growth is accounted for in one of three ways: either analytically, using 
dfilamination growth laws in conjunction with strain energy release rate analyses 
i ncorporat i ng delamination resistance curves; experimentally, using measured 
stiffness loss; or conservatively, assuming delamination onset corresponds to 
catastrophic delamination growth. Fail-safety is assessed by accounting for the 
accumulation of delaminations through the thickness. A tension fatigue life 
prediction for composite laminates is presented as a case study to illustrate 
how this approach may be implemented. Suggestions are made for applying the 
Damage -threshold/Fail-safety approach to compression fatigue, 
tension/compression fatigue, and compression strength following low velocity 
impact . 

Keywords: Damage Tolerance, Threshold, Fail-safe, Composite Materials, 
Delamination, Impact, Fatigue, Compression, Strain Energy Release 
Rate, Fracture Mechanics 
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NOMENCLATURE 


A Coefficient in power law for delamination growth 

a Delamination size 

b Laminate half-width 

c uncracked ply thickness 

d cracked ply thickness 

E Axial modulus of a laminate 

E . Axial modulus before delamination 
LAM 

* 

E Modulus of an edge delaminated laminate 

E Modulus of a locally delaminated cross section 
LD 

E Modulus of local cross section with edge and local delaminations 
L D 

E Initial modulus (N=10^ cycles) 
o 

E 1 1 Lamina modulus in the fiber direction 

E ^2 Lamina modulus transverse to the fiber direction 

G.^ In-plane shear modulus 

0 Strain energy release rate 
G Mode I strain energy release rate 

G Mode II strain energy release rate 

G Critical value of G at delamination onset 
c 

G Maximum G in fatigue cycle 
max 

K Strain concentration factor 
e 
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I Laminate length 

M Number of sublaminates formed by edge delamination 
m Slope of G versus log N curve for delamination onset 
n Exponent in power law for delamination growth 
N Number of fatigue cycles 
Np Cycles at failure in fatigue 

p Number of local delaminations through the laminate thickness 

R Cyclic stress ratio in fatigue (o . /a ) 

min max 

2s Matrix crack spacing 
t Thickness 

t, ... Laminate thickness 
LAM 

t Thickness of a locally delaminated cross section 
Ll) 

e Uniaxial strain 

e Critical strain at delamination onset 
c 

z Strain at failure 

r 

e Maximum strain in fatigue cycle 
max 

o Uniaxial stress 

o Maximum stress in fatigue cycle 
max 

o . Minimum stress in fatigue cycle 
min 

°alt A l ternatin 8 Stress in fatigue cycle 
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INTRODUCTION 


As composite materials are considered for primary structural applications, 
concern has been raised about their damage tolerance and long term durability. 
The threat of barely visible, low velocity impact damage, and its influence on 
compression strength, has surfaced as the most immediate concern for primary 
structural components such as composite wings [1]. Recent government programs 
have focused heavily on this issue in developing damage tolerance criteria that 
will satisfy the safety requirements of current military aircraft [2,3]. At the 
3ame time, research has been conducted on low velocity impact; both in the 
prediction of damage accumulation during the impact [H,5], and in the assessment 
of the influence of impact damage on compression strength [ 6 - 13 ]. Several 
methods for improving the performance of impacted composite panels and 
components have been proposed. One approach is to increase the inherent 
toughness of the composite by using tougher resin matrices, such as toughened 
epoxies [9] and thermoplastics [10], or to modify the form of the material by 
adding tough adhesive layers during the layup or as interleaves in the prepreg 
[12]. In terms of wing skin design, the goal has been to increase the 
compression failure strain after impact above the strength of a comparable 
laminate with an open hole [6,7]. Although this goal may be achieved using 
clever structural design and the improvements in materials cited, other issues 
have yet to be adequately addressed. 

Although compression strength is greatly reduced after low velocity impact, 
any further reduction with subsequent fatigue cycles is minimal. Hence, impacted 
composite panels have very flat compression S~N curves [1,6,13]. This 
observation has resulted in damage tolerance criteria for composite structures 
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that require only static loading [2]. However, for toughened matrix composites, 
where the compression strength after low velocity impact exceeds the strength of 
the laminate with an open hole, a static criteria may no longer be sufficient. 
The compression S-N curve for composite laminates with an open hole is not flat, 
even for toughened matrix composites [14], because the interlaminar stresses at 
the hole boundary cause delaminations that form in fatigue and grow with 
increased cycles [15]. Furthermore, other sources of delamination (straight 
edges, ply drops, matrix cracks) may exist in wing skins and other composite 
primary structures, such as composite rotor hubs [16]. Although delamination may 
not cause immediate failure of these composite parts, it often precipitates 
component repair or replacement, which inhibits fleet readiness, and results in 
increased life cycle costs. Furthermore, delaminations from several sources may 
accumulate, eventually leading to catastrophic fatigue failures. 

In metallic structures, damage tolerance has been demonstrated using 
fracture mechanics to characterize crack growth under cyclic loading for the 
constituent materials, predict the rate of crack growth in the structure under 
anticipated service loads, and establish inspection intervals and nondestructive 
test procedures to ensure fail safety. Because composite delamination is a 
commonly observed damage mechanism in laminated composite structures, many 
efforts have been undertaken to develop similar procedures for composite 
materials by characteri zing delamination growth using fracture mechanics [17— 
20]. Although this approach is promising, there are some fundamental differences 
in the way fracture mechanics characterization of delamination in composites may 
be used to demonstrate fail safety compared to the classical damage tolerance 
treatment used for metals. 

Previously, a Damage-threshold/Fai 1-saf ety approach to composite damage 
tolerance was proposed as an alternative to the classical approach used for 
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metals [21]. The purpose of the current paper is to expand on this concept by 
demonstrating how a Damage-threshold/Fail-safety approach may be used to predict 
the tension fatigue life of composite laminates, and then illustrating the 
similarities between this application and the use of the same philosophy for 
predicting compression fatigue life and compression strength after low velocity 
impact . 

DELAMINATION CHARACTERIZATION 


Many papers have been published recently where the rate of delamination 
growth with fatigue cycles, da/dN, has been expressed as a power law 
relationship in terms of the strain energy release rate, G, associated with 
delamination growth [17-20]. This fracture mechanics characterization of 
delamination growth in composites is analagous to that of fatigue crack growth 
in metallic structures, where the rate of crack growth with cycles is correlated 
with the stress intensity factor at the crack tip. However, delamination growth 
in composites occurs too rapidly over a small range of load, and hence G, to be 
incorporated into a classical damage tolerance analysis for fail safety 
[18,21,22], Where in metals the range of fatigue crack growth may be described 
over as much as two orders of magnitude in G, the growth rate for a delamination 
in a composite is often characterized over barely one order of magnitude in G. 
Hence small uncertainties in applied load may yield large (order of magnitude) 
uncertainties in delamination growth rates. 

Different damage mechanisms may also interact with the delamination and 
increase the resistance to delamination growth. Delamination growth resistance 
curves may be generated to characterize the retardation in delamination growth 
from other mechanisms [23-25]. These delamination resistance curves are 
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analogous to the R-curves generated for ductile metals that account for stable 
crack growth resulting from extensive plasticity at the crack tip. However, 
unlike crack tip plasticity, other composite damage mechanisms, such as fiber 
bridging and matrix cracking, do not always retard delamination growth to the 
same degree. Hence, the generic value of such a char act er i za t i on is 
questionable . 

One alternative to using the classical damage tolerance approach for 
composites as it is used for metals would be to use a strain energy release rate 
threshold for no delamination growth and design to levels below this threshold 
for infinite life. Metals are macroscopically homogeneous, and the initial 
3 tress singularities that create cracks at particular locations in preferred 
directions cannot be easily identified. Composites, however, are macroscopically 
heterogeneous, with stiffness discontinuities that give rise to stress 
singularities at known locations such as straight edges, internal ply drops, and 
orthogonal matrix cracks. Although these singularities are not the classical 
- 1 /p 

r variety observed at crack tips, and hence cannot be characterized with a 
single common stress intensity factor, they can be characterized in terms of the 
strain energy release rate, G, associated with the eventual delamination growth. 

The moat common technique for characterizing delamination onset in 
composite materials is to run cyclic tests on composite specimens, where G for 
delamination growth is known, at maximum load or strain levels below that 
required to create a delamination under monotonic loading. A strain energy 
release rate threshold curve for delamination onset may be developed by running 
tests at several maximum cyclic load levels and plotting the cycles to 
delamination onset versus the maximum cyclic G, corresponding to the maximum 
cyclic load or strain applied [26-30]. This G threshold curve may then be used 
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to predict delamination onset in other laminates of the same material, or from 
other sources in the same laminate [31]. 


DAMAGK -THRESHOLD/FAIL-SAFETY APPROACH 


One concern with a no-growth threshold design criteria for infinite life 
ha3 been the uncertainty inherent in predicting service loads, which could lead 
to G values that exceed no-growth thresholds and result in catastrophic 
propagation. This concern is paramount for military aircraft and rotorcraft, 
where original mission profiles used to establish design loads are often 
exceeded once the aircraft is placed in service. However, unlike crack growth in 
metals, catastrophic delamination growth does not necessarily equate to 
structural failure. In situations where the structure experiences predominantly 
tensile loads, such as composite rotor hubs and blades, delaminated composites 
may have inherent redundant load paths that prevent failure and provide a degree 
of fail safety [ 21 ]. This degree of fail safety has led some designers to think 
of composite delamination as a benign failure mode. Unfortunately, delaminations 
rnay occur from several sources in a given component or structure. When this 
occurs, an iterative composite mechanics analysis that considers each of these 
potential sites must be performed to ensure fail safety of the structure. 

Previously, a Damage-threshold/Fail-saf ety approach for composite fatigue 
analysis was proposed [21] that involved the following steps: 

(1) Predict delamination onset thresholds using fracture mechanics 

( 2 ) Assume complete propagation ocurrs immediately after the delamination 
threshold is exceeded 
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(3) Determine the remaining load carrying capability of the composite with 
delamination present using composite mechanics (i.e., check for fail 
safety) 

(4) Iterate on steps 1-3 to account for multiple sources of delamination 

This type of analysis need only be applied to primary structures. However, step 
1 may be used to demonstrate the delamination durability of any composite 
structure by providing an assessment of component repair or replacement costs 
over anticipated structural service lives. Step 2 reflects a conservative way to 
deal with the rapid delamination growth rates observed relative to metals as 
discussed earlier. An alternative to step 2 would be to predict delamination 
growth rates using growth laws that incorporate R-curve characterizations, 
thereby taking into account the resistance provided by other damage mechanisms. 
Such a characteri zation has been attempted previously [25], but should be used 
with caution because it is no longer truly generic. A third approach is to 
monitor stiffness loss in real time, and hence reflect the consequence of 
delamination growth, and other damage mechanisms, as they occur. This technique 
was used to predict the tension fatigue life of composite laminates [31 L and is 
summarized in the next section to provide a specific case study for the 
implementation of the Damage-threshold/Fail-saf ety approach. However, in most 
structural applications real-time monitoring of stiffness loss may not be 
practical, so the conservative approach outlined in step 2 would be applied. 
Step 3 acknowledges that the residual strength of the composite is a function of 
structural variables, and is not uniquely a question of material 
characterization. Hence, the Damage-threshold/Fail-saf ety concept offers both 
the benefits of generic material characterization using fracture mechanics, 
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while reflecting the unique structural character of laminated composite 


"materials 


LAMINATE FATIGUE LIFE PREDICTION: A CASE STUDY 


Tension Fatigue Behavior 


Figure 1 shows the tension fatigue damage in (45/-45/0/90) s X751/50 E-glass 

epoxy laminates that were subjected to cyclic loading at a frequency of 5 Hz and 
an R ratio of 0.1 [31]. Figure 2 shows a schematic of some of this damage, 
including edge delaminations that form at the edge in the 0/90 interface and 
jump through 90 degree ply cracks to the other 0/90 interface, and local 
delaminations that form in the 45/-45 interface, originating at 45 degree matrix 
ply cracks. These same damage mechanisms have been observed in graphite epoxy 
laminates with the same layup [27,35]. The fatigue damage in the glass epoxy 
laminates progressed in the following sequence as they were tested at maximum 
cyclic stress levels below their static strength. First, extensive matrix 
cracking developed in the 90 degree plies, followed by edge delamination in the 
0/90 interfaces. Next, matrix cracks appeared in the 45 degree and -45 degree 
plies and initiated local delaminations, first in the 45/-45 interfaces, 
followed by the -45/0 interfaces. Finally, after enough local delaminations had 
formed through the thickness at a particular location, fiber failure occurred 
and the laminate fractured. 

Fig. 3 shows the number of cycles at a given maximum cyclic stress for edge 
delaminations to form (solid symbols), for the first local delamination to form 
at the 45/~45 interface (brackets), and for fatigue failure to occur (open 
symbols). Under monotonic loading, matrix cracks formed in the 90 degree plies, 


10 



followed by edge delaminations in the 0/90 interfaces (solid symbols), and 
finally, by fiber fracture (open symbols). In order to predict the ultimate 
fatigue failure of these laminates, the onset and growth of the damage observed 
must be characterized, and the influence of this damage on laminate stiffness 
and strength must be determined. Once these relationships are known, fatigue 
life may be predicted using the the Damage-threshold/Fail-saf ety approach. 

Influence of Damage on Laminate Stiffness 


Fig. 4 shows the influence of damage on laminate stiffness. As matrix 
cracks accumulate, and as delaminations form and grow, the stiffness of the 
laminate decreases. Laminate stiffness is the ratio of the remote stress to the 
global strain in the laminate. The global strain is typically measured using an 
extensometer or LVDT (fig. 4), which yields the displacement of the laminate over 
a fairly long gage length relative to the laminate’s length. As damage forms and 
grows in the laminate under a constant maximum cyclic stress, corresponding to a 
constant applied maximum cyclic load, the global strain in the laminate 
increases . 

Previous studies have determined the relationships between stiffness loss 
and damage extent [23. 32-35 ]. The amount of stiffness loss associated with 
matrix cracking depends upon the ply orientation of the cracked ply, the 
laminate layup, the relative moduli of the fiber and the matrix, and the crack 
spacing, or density of cracks, in the ply. For example, in ref. 33. an equation 
was derived for stiffness loss due to matrix cracking in the 90 degree plies of 
cross ply laminates as 
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E 


( 1 ) 


e lam 

1 + (1/Xs)(c/d)(E 22 /E 11 )tanh(As) 


where 


X 


( 3G 12 (G+d)E LAM j 
c 2 dE u E 22 


1/2 


( 2 ) 


and c and d are the thicknesses of the cracked and uncracked plies, 
respectively. As the crack density increases, i.e., as the crack spacing, 2s, 
decreases, the stiffness of the laminate will decrease. 

The amount of stiffness loss due to delamination also depends on the 
laminate layup and the relative moduli of the fiber and the matrix, as well as 
the location and extent of the delamination. As delaminations form and grow in a 
particular interface, the laminate stiffness decreases as the delamination size, 
a, increases. In ref. 23, an equation was derived for the stiffness loss 
associated with edge delamination as 

E ■ <E ‘- E LA« )a/b * E LAM <3) 

where a/b is the ratio of the delamination size to the laminate half-width, and 
* 

E is determined from a rule of mixtures expression 


* 

E 


M 

I 

i = 1 


E i t i 


(4) 
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are 


where the moduli of the M sublaminates formed by the delamination, 

* 

calculated from laminated plate theory. The difference in E, ... and E reflects 

LAM 

the loss of transverse constraint in the sublaminates formed by the 
delamination . 

Delaminations starting from matrix cracks will affect laminate stiffness 
differently than delaminations growing from the straight edge. In ref. 32, an 
equation was derived for the stiffness loss associated with delaminations from 
matrix cracks as 


= ( (a/JOt, . M (1/t. E - 1/t tAM E, AM ) 
LAM LD LD LAM LAM 


»-1 


( 5 ) 


where a/l is the ratio of the delamination length to the laminate length, and 

E and t represent the modulus and thickness of the locally delaminated 
L D LD 

region in the vicinity of the matrix crack. The locally delaminated modulus, E. , 

LD 

* 

is calculated using laminated plate theory and is similar to E in eqC 1 !). 
However, in addition to reflecting the loss in transverse constraint due to the 
delamination, E LQ also reflects the loss of the load bearing capacity of the 

cracked ply. Similar to edge delamination, the stiffness of the laminate 
decreases as the size of the delamination increases. However, unlike edge 
de Laminations , which form at the two edges and grow progressi vely towards the 
center of the laminate width, local delaminations tend to accumulate at several 
matrix cracks along the length, growing only a small distance at any one 
location. The cumulative effect of these local delaminations with cycles, 
however, may have a significant effect on measured stiffness loss. 
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Delamination Onset and Growth Characterization 


In order to predict stiffness loss as a function of fatigue cycles, the 
onset and growth of matrix cracks and delaminations must be characterized. This 
characterization should be done using a generic parameter that is representative 
of the composite material being tested, but independent of laminate structural 
variables such as layup, stacking sequence, and ply thickness. Typically, 
the strain energy release rate, G, associated with matrix cracking and 
delamination is used for this purpose. 

Fig. 5 shows the steps that would be required to predict stiffness loss as 
a function of fatigue cycles using a G characterization of damage onset and 
growth. First, plots of G versus log N must be generated to characterize the 
onset of delamination [21,28,31], and power law relationships between G and the 
rate of growth of delamination with fatigue cycles are needed to characterize 
damage growth [ 17 - 20 , 22 ] (fig. 5a). Similar approaches may be used to 
characterize the onset and growth of matrix cracks [33]- Using these material 
characterizations, the decrease in matrix crack spacing, 2s, and the increase in 
delamination size, a, with fatigue cycles may be predicted (fig. 5b). This 
information, in turn, may be used with eqs(1-5) to predict the decrease in 
modulus with cycles, which for a constant stress amplitude test is tantamount to 
predicting the increase in global strain with cycles (fig. 5c). 

Although this technique may be demonstrated for cases where there is one 
dominant damage mechanism, application of this approach in general is difficult 
because the various damage modes interact, complicating their unique 
characterization in terms of G [33]. For example, although the elastic analysis 
for G associated with edge delamination growth is independent of delamination 
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size, stable delamination growth is often observed experimentally [23,25]. The 
strain energy release rate for edge delamination was derived in ref. 23 as 


G 


A 


LAM 


(e lam- e *> 


( 6 ) 


which is independent of the delamination size. Theoretically, when a critical 

value of strain, e , is reached, corresponding to a critical G , the 

c c 

delamination should grow catastrophically across the laminate width. However, 

edge delaminations usually grow in a stable fashion, requiring increasing strain 

levels, and hence increasing G, for the delamination to grow across the width. 

This stable growth may be correlated with the accumulation of 90 degree matrix 

cracks ahead of the delamination front. For example, fig. 6 shows a plot of 

normalized delamination size, a/b, as a function of the strain applied to an 

eleven-ply ( ±30/±30/90/90) T300/5208 graphite epoxy laminate [23]. The edge 

3 

delaminations form at a strain of approximately 0.0035, but do not grow across 
the specimen width until the strain reaches approximately 0.0065. Also plotted 
in fig. 6 on the right hand ordinate is the 90 degree ply crack spacing measured 
in the center of the laminate. There appears to be a direct correlation between 
the stable delamination growth and the accumulation of matrix cracks ahead of 
the delamination front. These matrix cracks apparently alter the local stress 
state at the delamination front and increase the resistance to delamination 
growth. Stable delamination growth may be predicted by generating a delamination 
resistance R-curve using eq.6 [23*25]. However, the resulting R-curve is no 
longer generic, because the matrix cracking that is causing the delamination 
resistance is governed by structural variables such as ply thickness and 
stacking sequence. 
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Delamination also influences the formation and accumulation of matrix 


cracks. Delamination relaxes the constraint of neighboring plies, and hence 
changes the saturation spacing of matrix cracks in the off-axis plies. For 
example, when delaminations form at the edges of the (±30/±30/90/90) laminate 

3 

3hown in fig. 6, the constraint between the -30 deg and 90 deg plies is relaxed, 
and the 90 deg cracks form at lower strain levels, with smaller crack spacings, 
than possible if no delamination had existed [23]. An R-curve description of 
matrix cracking has been used to describe the accumulation of matrix cracks, 
similar to the approach that has been attempted for delamination [331- However, 
when these cracks interact with delaminations, this description is no longer 
generic . 

Even if one could achieve a truly generic description of damage 
accumulation with cycles, the resulting stiffness loss prediction, and hence the 
prediction of increasing global strain with cycles, is necessary, but not 
sufficient, to predict fatigue life. The final failure of the laminate is 
governed not only by loss in stiffness, but also by the the local strain 
concentrations that develop in the primary load bearing plies, which in most 
laminates are zero degree plies. 

Influence of Local Strain Concentrations on Failure 


Fig. 7a shows that fatigue failures typically occur after the global strain 
has increased because of tne fatigue damage growth, but before this global 
strain reaches the global strain at failure, , measured during a static 

strength test [27,31,3^,35]. Therefore, local strain concentrations must be 
present in the zero degree plies that control the laminate strength. Although 
matrix cracks create small strain concentrations in the neighboring plies, their 
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magnitudes are generally small because the stiffness of the cracked ply is 
usually much less than the stiffness of the zero degree ply [36]. Furthermore, 
strain concentrations due to matrix cracks act over only a local volume in the 
adjacent ply near the crack tip [ 36 ]. Hence, the final failure in a zero degree 
ply of a laminate may follow a neighboring ply crack [37], but the laminate 
failure strain will not be strongly influenced by the presence of the matrix 
cracking [ 36 ]. However, once delaminations initiate at matrix ply cracks 
anywhere through the laminate thickness, the local strain will increase 
significantly in all of the remaining uncracked plies [21 ,27, 32,3^ — 363. These 
local strain increases may not be sensed by the global strain measurement, 
because delaminations starting from matrix cracks grow very little once they 
form. If several delaminations form at matrix cracks throughout the laminate 
thickness at one location, then the local strain on the zero degree plies at 
that location may reach the static failure strain, resulting in the observed 
fatigue failure (fig. 7b). 

This mechanism for fatigue failure has been observed previously for 
graphite epoxy laminates [27,3^.35], and the local strain concentrations 
resulting from cumulative local delaminations through the thickness have been 
quantified [21 , 27 , 32 , 3 1 * ■ 35] . These local strain concentrations may be calculated 
as 


K 

z 


E LAM t ’LAM 

E LD t LD 


( 7 ) 


Typically, the local strain concentration will result in a trade off between the 

increased modulus, E, E, because E, ^ is a more zero degree dominated layup 

LD LAM LD 
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than the original laminate, and the decrease in load bearing cross section, 
t LD < t LAM - 

Each time a delamination initiates from a matrix crack, the local strain in 

the remaining uncracked plies, and hence in the zero degree plies, increases by 

an amount equal to K times the global cyclic strain, e , until it reaches the 

e max 

static failure strain, e , (fig. 8a). A simpler way to visualize this process, 

F 

however, is to reduce the static failure strain to some effective global failure 
strain, (£„)., each time a new local delamination forms through the thickness. 

r 1 

Hence, (c„). would be equal to e„/(K ).. As local delaminations accumulate 
F 1 F e x 

through the thickness, the effective failure strain would decrease 
incrementally. Because (K ^ reflects the strain concentration resulting from 

accumulated local delaminations through the laminate thickness, fatigue failure 

would correspond to the number of cycles where the damage growth increased the 

global maximum cyclic strain, e , to the current value of (e„). (fig. 8b). 
° max F x 

This approach does not require a prediction of damage growth with fatigue cycles 

if the laminate stiffness los3, and hence the increase in global strain, can be 

monitored in real time. When this is possible, only the incremental decrease in 

the effective failure strain needs to be predicted to predict fatigue life. This 

may be accomplished by assuming that matrix cracks exist in all of the off-axis 

plies. This assumption is analagous to assuming the existence of the smallest 

flaw in a metal that could be detected non-destruct i vely to assess damage 

tolerance. Then, the number of fatigue cycles to onset of each local 

delamination through the thickness may be predicted using delamination onset 

criteria (fig. 5a) along with strain energy release rate analyses for local 

del am i na t i on . As each local delamination forms, may be reduced by the 

r 
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appropriate K and compared to the current value of e , based on measured 
£ max 

stiffness I 033 , to determine if fatigue failure will occur. Hence, the ability 

to predict local delamination onset, and its effect on e , facilitates using 

F 

measured stiffness loss to predict fatigue life. However, for many composite 
structures real-time stiffness measurement may not be practical. In these cases, 
the conservative approach for step 2 in the Damage-threshold/ Fail-safety 
approach outlined earlier could be applied. 

If the conservative approach was used to predict the tension fatigue life 
of ( 45/ — ^5/0/90) laminates, for example, stiffness would decrease 

3 

incrementally, i.e., e would increase incrementally, with the onset of each 

ms x 

damage mechanism. Figure 9 shows a sketch for conservative fatigue life 

prediction in ( 45/-45/0/90) graphite-epoxy and glass-epoxy laminates. Because 

s 

matrix cracks are assumed to exist in the off-axis plies, e is increased in 

K max 

the first load cycle, corresponding to the stiffness loss associated with 
saturation crack spacing in the off-axis plies. This stiffness loss would be 
greater for glass epoxy laminates than for graphite epoxy laminates [30. The 
influence of matrix cracks on local strains in the zero deg plies will be 
neglected for the reasons stated earlier. Hence, e^. will remain unchanged. When 

edge delamination occurs in the 0/90 interfaces, e will increase again, 

cor r espond i ng to complete delamination throughout the laminate width. This 
stiffness loss would be greater for graphite epoxy laminates than for glass 
epoxy laminates [30. However, e would not change because edge delaminations do 

r 

not create local strain concentrations in the zero deg plies [32]. As each local 

delamination forms, (e ). will decrease, as determined by (K )., and c will 

F l e l max 

increase, corresponding to delamination growth throughout the particular 
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interface. When enough local delaminat ions form through the thickness such that 

(e ). 2 (e„). t fatigue failure will occur. These predictions will be 
max 1 r 1 

conservative because matrix cracking typically does not reach saturation spacing 

in all of the off-axis plies in the first load cycle, and because 

delaminations grow rapidly at first, but then are retarded by interaction with 

matrix cracking, and hence, rarely grow catastrophically. 

Because stiffness loss data were available for the glass epoxy laminates in 

this case study, measured stiffness loss was used to determine the increase in 

c with fatigue cycles instead of using the conservative prediction 
m 3 x 

methodology. Furthermore, the G vs. log N delamination characterization was 

generated using edge delamination data from the ( H5/-H5/0/9Q) laminates and was 

s 

then used to predict local delamination onset in these same laminates. Hence, 
the accuracy of this fatigue life prediction depends primarily on the validity 
of reducing the failure strain incrementally to account for the accumulation of 
local delaminations through the laminate thickness. The next section outlines 
how this fatigue life prediction was performed in the context of the Damage- 
threshold/ Fail-safety philosophy. 

Life prediction using Damage-threshold/Fail-safety approach 
Step 1: Delamination onset prediction 


In order to predict the onset of local delaminations with fatigue cycles, 
the G versus log N characterization of the composite material must be generated. 
This characterization may be accomplished using a variety of interlaminar 
fracture test methods [22,26,28-30]. Data from several materials with brittle 
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and tough matrices indicated that between 10^5 N £ 10^ cycles, the maximum 

cyclic G may be represented as a linear function of log N (fig. 10), where N is 

the number of cycles to delamination onset at a prescribed G [28]. Hence, 

max 


G=mlogN+G c (8) 

where G and m are material parameters that characterize the onset of 

delamination under static and cyclic loading in the material (fig. 5a). A more 
recent study indicates that for some tough thermoplastic matrix composites, the 
static toughness is so great compared to the fatigue behavior that a linear 
representation may not be valid [22]. Furthermore, the static toughness may 
depend on the loading rate, which would influence the ultimate shape of the 
curve. For example, if slow monotonic loading is used to measure the static 
toughness, both G c and the slope, m, may be different from the results plotted 

using G c measured at a load rate that corresponds to the frequency of cyclic 

loading. Because a low toughness, glass epoxy composite was tested in ref. 31, a 
linear characteri zation was used in this case study. The static data in ref. 31 
were measured at a relatively slow rate of 0.5 mm/min. 

To predict delamination onset, G must be calculated for the first local 
delamination that will form. This typically occurs at a matrix crack in the 
surface ply, but may be confirmed by calculating G for matrix cracking in all of 
the off-axis plies in the laminate. The one with the highest G for the same 
applied load will be the first to form. This G may be calculated using the 
equation for the strain energy release rate associated with local delaminations 
initiating at matrix cracks [32] 


21 



G 


(9) 


2 2 

0 ^AM 


< 1 /t LD^'LD - 1/t LAM E LAM ) 


To calculate the number of cycles for the first local delamination to form, 

eq(9) for G is set equal to the delamination onset criterion of eq(8) and then 
solved for . Hence, 


log N, - - 


[ JMX t^ AM ( 


1/t, „E, - 1/t, AM E r 1 - G ] 
LD LD LAM LAM'' C J 


( 10 ) 


Step 2: Assessment of Damage Growth and Stiffness Loss 


Delamination growth information is needed to determine the amount of 
stiffness 1033 , and hence the increase in global strain, that ha3 occurred by 
the time the first local delamination has formed at N 1 cycles. In graphite epoxy 

laminates, the majority of this stiffness loss is associated with delamination; 
however, in glass epoxy laminates matrix cracking may also contribute 
significantly to stiffness loss [31]. In either material the interaction of 
matrix cracking and delamination complicates the prediction of damage growth, 
and hence the prediction of stiffness loss. Therefore, instead of predicting 
stiffness loss by predicting the rate of delamination growth and accumulation of 
matrix cracks with fatigue cycles, stiffness loss was monitored experimentally. 

Step 3- Assessment of fail-saf ety 
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The strain concentration associated with the first local delamination, 
(K c ) lf may be calculated using eq.7. Fatigue failure will occur if the maximum 

global strain, resulting from the stiffness loss associated with damage growth 
at N cycles, reaches the effective failure strain when the local delamination 

forms, which is calculated as (e p ) 1 = e p /(K £ ) 1 . Hence, failure will occur if 

e „„ - (e-.).* If the first local delamination does not cause failure, then 
max r i 

further local delamination sites must be considered. 

Step 4: Analysis of multiple local delaminations through the thickness 


As shown in fig.11, the thickness and modulus terms in eq(9) change for 
each succesive local delamination that forms through the thickness. For example, 
t and E for a 45/-45 local delamination in a (H5/-H5/0) laminate becomes 

LL) LI) S 

the t LAM and E LAM value3 U3ed for ttie next l° ca I delamination that forms through 

the thickness. Therefore, as local delaminations accumulate through the 

thickness under a constant o , the driving force (i.e. G) for each new 

m 3 x 

delamination changes. Hence, fatigue life prediction for composite laminates 
requires a "cumulative damage" calculation, even for constant amplitude loading. 
To calculate the number of cycles for each successive local delamination to 
form, , the appropriate form of eq(9) for G is set equal to the delamination 

onset criterion of eq(8) (fig.11) and then solved for N. (fig. 12). Hence, 


log N 


i 


i [ 


max 






1/t, . m E. - G ] 

LAM LAM'l C J 
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Fatigue failure will occur when e Z (e„)., with a resulting fatigue life, N„, 

max F 1 F 


of 


P 

N = I N (12) 

i = 1 

where p is the number of local delaminations that form through the thickness of 
the laminate before failure. 

Because of the scatter in the experimental data, the constant load 

amplitude fatigue life prediction methodology outlined in fig. 12 more closely 

resembles fig. 13- The variation in initial laminate modulus (i.e, the variation 

in c ) and the variations in the static failure strains from specimen to 
max 

specimen must be taken into account. Hence, a range of possible fatigue lives 

would be predicted, rather than a single value. The lowest life would occur when 

the minimum value of (e„). in the e 0 distribution reaches the largest e value 

F i F max 

in the distribution resulting from variations in laminate moduli. 

Life Prediction for Glass-Epoxy Laminates 


The Damage-threshold/Fai 1-saf ety approach outlined above was used to 
predict the fatigue life of (45/-4 5/0/90 ) E-glass epoxy laminates [31]. First, 

the delamination onset behavior in fatigue was characterized in terms of strain 
energy release rates. The maximum cyclic strain versus cycles to edge 
delamination onset for the laminate was used in eq(6) and the data were plotted 
versus log N (fig. 14). There was significant scatter in the static data for G , 

V 

possibly due to the interaction that occurred between the edge delamination as 
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it formed and the 90 deg ply cracks that were extensive before edge delamination 

onset [33]- Previous work has demonstrated that G values from edge delamination 

c 

data may be artificially elevated if extensive 90 degree cracking is present in 

the laminate [26J. Therefore, the minimum values in fatigue were used in eq(8) 

to characterize delamination onset. For the X751/50 E-glass epoxy, a G value of 

c 

2 

0.56 in-lbs/in was obtained, and the slope, m, was -0.06. 

Figure 15 shows the maximum cyclic strain as a function of fatigue cycles 
for the ( 45/- 1 l5/0/90) s X751/50 E-glass epoxy laminates cycled at a maximum 

cyclic stress of 210 MPa and an R of 0.1. Also shown in fig. 15 is the reduction 
in effective e for local delaminations accumulating through the thickness. The 

r 

range of estimated and measured fatigue lives for several 0 levels is 

max 

summarized in fig. 16. The agreement between predicted and measured fatigue lives 
is reasonably good. 

FACTORS AFFECTING DELAMINATION ONSET AND GROWTH 


The agreement between measured and predicted fatigue lives in fig. 16 
indicates that the Damage-threshold/Fai 1-safety approach, in the form of a 
through-thickness damage accumulation model, can accurately describe fatigue 
failure for a material whose delamination behavior in fatigue is well 
character i zed . In this case, the G vs. log N characterization was generated 
using data from the same laminates whose fatigue lives were being predicted. In 
general, however, the G vs. log N characterization would be performed on 
standardized laboratory tests, and then used to predict the fatigue behavior of 
structural components made of the same material. Hence, the laboratory 
characterization must be performed on identical materials (same constituents, 
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fiber volume fraction, cure conditions, etc.) under identical environments 
(temperature, moisture, etc.) and loading conditions (load rate, R-ratio, 
frequency, etc.) as the structure for the fatigue life prediction to be 
accurate. Furthermore, although delamination growth data are difficult to 
utilize because of steep growth rates and damage mode interactions, these data 
are useful, nevertheless, to identify how the various material, environment, and 
loading variables that effect delamination onset will influence delamination 
growth . 

Of the many factors that may affect delamination onset and growth, a few 
have been studied in detail. For example, the toughness of the matrix will have 

a very stong effect on G c but very little influence on delamination onset at 10^ 

cycles (fig. 10) [26,28-30]. Therefore, the slope, m, as measured by fitting the 
delamination onset data to eq(8) will be lower for a brittle matrix composite 
than a tougher matrix composite (f ig. 1 7a)- [28] . Assuming that the brittle and 
tough matrix composites eventually reach a common G threshold for delamination 

onset at N£10^ cycles [22,26,28,29], then the exponent, n, in a delamination 
growth law would be lower for the tougher matrix composite (fig. 17b) [18]. 

Delamination characterization may also depend on the mixed-mode ratio for 
the particular source of delamination. Previous studies have shown that the 
total G c at delarnination onset under a monotonic loading varies as a function of 

the mixed-mode percentage at the delamination front [29,38,39]. The total G q 

will be highest for situations where the mode II component is greater than the 
mode I component (fig.l8a). However, the G threshold for delamination onset at 

N£ 1 0^ cycles ha3 been shown to be nearly identical for all mixed mode ratios, 
from pure mode I to pure mode II [22,26,28,29]. Therefore, as shown in fig.l8a, 
m in eq(8) will be greater for delaminations that are predominantly due to 
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interlaminar shear (mode II) than for delaminations that are predominatly due to 
interlaminar tension (mode I). Assuming a common G threshold for delamination 

onset at N£10^ cycles, the exponent in the delamination growth power law would 
be lowest for the pure mode II case and highest for the pure mode I case 
(fig.l8b). Previous delamination growth studies have verified these trends [IT- 
19]. For the glass epoxy laminate fatigue life prediction summarized earlier, 
conservative values of G c and m were used in eq(8) because of the scatter in the 

static total G c measured using edge delamination data. Hence, the mixed-mode 

ratio dependence was ignored. In general, however, the mixed mode dependence on 

G should be determined for both the material characterization test(s) as well 
c 

as the delamination source being modeled in the structural component. However, 
if the long term delamination durability is of primary concern, the G threshold 

at N£10^ cycles is all that is needed. In this case, only a simple total G 
analysis is required, since the G threshold does not depend strongly on the 
mixed mode ratio. This greatly simplifies the analysis, because total G may be 
calculated using relatively simple analyses like eqs(6) and (9) [23,32]. 

Changing the R-ratio of the cyclic loading will not affect G but may have 

a significant influence at 10^ cycles (fig. 19a) [22,30]. Therefore, the slope, 
m, will be greater for lower R-ratios corresponding to greater alternating 

stress levels. Hence, G threshold values at 10^ cycles will be lower for smaller 
R-ratios [30]. Consequently, the exponent of the delamination growth power law 
will be lower for the lower R-ratios (fig. 19b). 

The influence of other material, environmental, and loading variables have 
been examined [2*1, 40, Ml ] . However, most of this work has been performed for 
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static toughness and/or delamination growth. Much work still needs to be done to 
determine the influence of these variables on delamination onset. 


DAMAGE -THRESHOLD/FAIL-SAFETY APPROACH FOR COMPRESSION 


In the previous case study, and in the examples cited in ref. 21, the 
Damage-threshold/Fail-saf ety approach was illustrated for problems that involved 
only tension loading. However, this same approach may be applied to laminates 
subjected to compression loading. Delamination onset characterization would be 
conducted in the same way, with only the assessment of fail safety (step 3) 
changing significantly. 

The significance of accumulated delaminations on compression strength has 
been documented previously by comparing the strength of laminates with one, two, 
or three implanted delaminations through the thickness to identical laminates 
with either barely visible or visible impact damage (fig. 20) [42]. These results 
show that the compression strength for laminates with 2.U inch diameter 
implanted delaminations, normalized by the compression strength for the same 
laminates with a 1/4 inch open hole, decreases as the number of delaminations 
Increases through the thickness. Still lower compression strengths were observed 
for the impacted laminates, which typically contain delaminations in nearly 
every interface [11]. Similar studies have compared the residual compression 
strength of virgin laminates, or laminates that had implanted delaminations in a 
single interface, to identical laminates without implants that had undergone low 
velocity impact with subsequent cycling [6,13]. For example, fig. 21 shows a plot 
of cycles to failure as a function of stress amplitude for 
(0/90/0/45/-45/0) a graphite epoxy laminates subjected to fully reversed cyclic 

loading, either in the initially undamaged state, or following an impact with a 
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potential energy per unit thickness of 1790 J/m [ 13 ]. The data in fig. 21 
indicate that the compression strength after impact is very low compared to the 
fatigue behavior of the virgin laminate. Furthermore, most of the strength 
reduction occurs after the impact, with very little degradation due to 
subsequent cyclic loading. 

For composites loaded in compression, final failure is not necessarily 

determined by the local strain concentration in the zero degree plies, but often 

results from a global instability that occurs after delaminations accumulate 

through the thickness and become locally unstable. For example, fig. 22 shows 

dye penetrant enhanced radiographs of the edge of a MO-ply thick, 

( H5/0/-M5/90).. T300/3501-6 graphite epoxy laminate, containing Kevlar stitches 

5s 

across the specimen width, that was cycled in compression at a maximum cyclic 
compression stress of 32.5 Ksi and an R-ratio of 10 [433. After 320,000 cycles, 
delaminations had formed at the edge near the top surface. The sublaminate that 
formed buckled locally, which in turn led to more delaminations forming in 
adjacent interfaces and subsequently buckling. The accumulation of these 
delaminations through the thickness eventually reduced the cross section 
carrying the compression load to the point at which global instability occurred 
and the laminate fractured. This accumulation of del am i nat i ons through the 
thickness occurred over the last 1000 cycles of the fatigue life. In laminates 
without through-thickness stitching, this final phase of the fatigue life may be 
even more rapid, and very difficult to document. In these situations, where the 
accumulation of delamination through the thickness occurs rapidly, prediction of 
the initial delamination onset may provide a reasonable estimate of fatigue life 
in addition to establishing the delamination durability of the composite. 

Because of this progressive buckling mode of failure, compression fatigue 
lives are typically much lower than tension fatigue lives for identical 
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laminates subjected to identical load amplitudes [13]. Combined 
tension/compression fatigue lives may be reduced even further as a result of 
delaminations forming from matrix cracks under tension loads and then growing as 
a result of local instabilities under the compression loads [HH], In each case, 
however, the final failure results from an accumulation of delaminations through 
the thickness. The Damage-threshold/Fail-saf ety approach could be used to 
estimate fatigue lives in each case. First, delamination onset would be 
predicted using the appropriate analysis for G in eq.(8) depending upon the 
source of the original delamination. Next, delaminations would be assumed to 
grow throughout the interface immediately, or solutions for instability driven 
delamination growth in compression would have to be incorporated if stiffness 
loss could not be monitored directly in real time. Several fracture mechanics 
models have been developed for the growth of through-width and elliptical patch 
delaminations in a single interface [ 1 »5“ 1 »9]. These analyses would have to be 
extended to model laminates with multiple edge delaminations to simulate 
compression fatigue damage, and laminates with multiple delaminations that were 
formed by matrix cracks to simulate tension/ compression fatigue damage. Finally, 
fail safety may be assessed in compression, as delaminations form near the 
surface and then accumulate through the thickness, using appropriate models for 
local and global buckling of the damaged laminate. 

These same models could be used to evaluate the consequence of low velocity 
impact damage. Previous studies have shown that low velocity impact damage 
develops as extensive matrix cracking and associated delaminations through the 
thickness [10-12]. Delamination onset in these cases has been modeled as 
delaminations initiating from matrix cracks under bending loads [50]. In brittle 
matrix composites, impacts that are barely visible on the impacted surface may 
be extensive not only on the back surface, but throughout the laminate 
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thickness. This extensive deiamination results in greatly reduced compression 
strength. Subsequent cyclic loading may create only slightly greater damage 
growth, which would explain the relatively flat S-N curves observed for impacted 
brittle matrix laminates (fig. 21). Tougher matrix composites, however, suppress 
some of the delaminations that would otherwise form through the thickness during 
the impact [10]. Therefore, the compression strength following impact is greater 
than the compression strength for similar laminates with brittle matrices, but 
subsequent cyclic loading may cause further damage and corresponding reductions 
in residual compression strength. In either case, the Damage-threshold/Fai 1- 
safety approach may be used to characterize the delamination onset and assess 
the fail safety of the damaged laminate. 

SUMMARY 


© A Damage--threshold/Fai 1-safety approach was proposed to ensure that 
composite structures are both sufficiently durable for economy of 
operation, as well as adequately fail safe or damage tolerant for flight 
safety. This approach invloved the following steps: 

1) Matrix cracks are assummed to exist throughout the off-axis plies 

2) Delamination onset is predicted using a strain energy release rate 
char act eri zation 

3) Delamination growth is accounted for in one of three ways: 

0 Analytically, using delamination growth laws in conjunction with 
strain energy release rate analyses incorporating delamination 
resistance curves 

0 Experimentally, using measured stiffness loss 
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© Conservatively, assuming delamination onset corresponds to 
catastrophic delamination growth. 

4) Fail-safety is assessed by accounting for the accumulation of 
delaminations through the thickness. 

© A tension fatigue life prediction for glass epoxy laminates was 
presented a3 a case study to illustrate how the Damage-threshold/Fai 1- 
safety approach may be implemented. A fracture mechanics analysis of 
edge delamination was used to generate a delamination onset criterion 
for the material. Then, strain energy release rates were calculated for 
local delaminations that formed at matrix ply cracks through the 
laminate thickness, and where compared to the criterion to predict local 
delamination onset. Delamination growth was accounted for experimentally 
using measured stiffness loss. Finally, fail safety was determined by 
accounting for the local strain concentration on the zero degree plies 
resulting from delaminations forming at matrix cracks through the 
laminate thickness. 

© Suggestions were made for applying the Damagc-threshold/Fail-safety 
approach to compression fatigue, tension/compression fatigue, and 
compression strength following low velocity impact. In all of these 
analyses, strain energy release rates may be used to predict 
delamination onset, and fail safety may be assessed by accounting for 
the effect of delaminations that have accumulated through the thickness. 
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Figure 1. - Glass -epoxy tensile fatigue test. 
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Figure 2. - Schematic of delamination in quasi -isotropic laminate subjected 
to cyclic tension. 



43 


Figure 3. - Tension fatigue behavior of (45/-45/0/90) s X751/50 E-glass 
epoxy laminates. 
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Influence of damage on laminate stiffness. 
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Figure 5. - Prediction of stiffness loss in composite laminates. 
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accumulation of 90-degree ply matrix cracks. 


strain at failure in static test 
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Figure 7. - Increase in global versus local strain in zero degree plies. 
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(a) Graphite epoxy (b) Glass epoxy 

Figure 9. - Damage-threshold/fai 1 -safety analysis. 
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Figure 10. - Mechanical strain energy release rate at del ami nation onset as 
a function of fatigue cycles. 
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Figure 11. - Strain energy release rate for local del ami nation onset. 
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Figure 12. - Reduction in global ep due to local delamination 
accumulation through the thickness. 
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Figure 15. - Tension fatigue life estimation for (45/-45/0/90) s X751/50 
E-glass epoxy. 



[±45/0/90] X751/50 E-glass Epoxy 



56 


Figure 16. - Fatigue life prediction based on local delamination 
accumulation through the thickness. 






o 

57 


Figure 17. - Effect of matrix toughness on delamination onset and growth. 
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Figure 18. - Effect of mixed-mode ratio on delamination onset and growth. 
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Figure 19. - Effect of R-ratio and mean load on delamination criteria. 
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Cycles to failure 

Fiqure 21. - Fatigue behavior of initially undamaged and impacted graphite 
epoxy laminates under fully reversed cyclic loading. 
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Figure 22. - Radiograph of through-the-thickness damage. 
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